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{. INTRODUCTION 


A. PROBLEM DESCRIPTION 

The United States must be able to deploy U.S. Forces and material in support of 
a Commander conducting operations anywhere in the world. It is imperative that the 
forces and material arrive at the rime and location requested. Furthermore, for a 


deplovment to succeed the air and Sealift assets must be properly allocated. 


B. JOINT DEPLOYMENT AGENCY 

The complexity and magnitude of deploving U.S. Forces to an overseas area 
requires careful and thorough coordination. In March 1979 the Joint Deployment 
Agency (JDA) was established to be the single point of contact for strategic 
deployment planning and coordination. 

The Joint Deployment Agency’s mission 1s to support the Joint Chiefs of Staff 
(JCS) and Supported Commanders in planning for and executing deplovments. As 
directed by the Joint Chiefs of Staff, the JDA is responsible for coordinating the 
actions of deploying units and the land, air and sealift assets to be used in a joint 
deployment. The JDA also serves as the focal point for information associated with 
deployment decisions. 

The land, air and sealift assets to be used in a joint deployment are controlled by 


three different organizations which are called the transportation operating agencies. 


The Military Sealift Command (MSC), which operates some 150 ships in support 
of the Navy’s daily operations, has earmarked 51 additional ships for joint 
deplovment. These are comprised of nine dry cargo ships, 22 general-purpose oil 
tankers, 12 maritime prepositioning ships now on station at various points 
around the globe and eight fast sealift ships. In wartime, this number would be 
Supplemented on short notice by | lo ships in the Ready Reserve Fleet (RRF). 


The Mihtary Airlift Command (MAC), which is responsible for all of the 
mulitary’s overseas air and freight transportation, will immediately contribute 234 
C-141 “Starlifters” and 77 C-5 aircraft. In a protracted conflict, MAC will also 
Meanerceesse1O, 256 aircraiu 1m the Civil Reserve Air Fleet. Not included in these 
totals are some 500 of MAC’s C-130 “Hercules” transports that are meant for 
supporting overseas unified commands. 


The Military Traffic Management Command (MTMC) 1s responsible for 
idenufving and managing the transportation routes to be used in the United 
States for the shipment of supplies to the ports where they are to be loaded 
aboard MSC’s ships. MTMC will supervise the loading operations in the United 
States and the unloading at the destination. [Ref. 1: p. 40] 


C, DEPLOYMENT PEANRING 
1. Deliberate Planning 

Deplovment planning is normally a time-consuming and “deliberate” process. 
During an emergency or time-sensitive situation the JDA must be prepared to conduct 
deplovment planning in a “crisis” mode. The deliberate planning process takes place 
over a period of weeks or months. The procedure insures people and material can be 
moved to support military objectives. Every movement plan goes through a refinement 
process. The most important question in the process is, “Will the deployment plan 
work?” Once it is determined the deployment plan can be supported with available 
transportation assets, the JDA looks for ways to improve the plan. Can the 
deployment be carried out faster, more efficiently, more cost-effectively, or withmless 
vulnerability to enemy attack? Refinement is an ongoing process with movement 
requirements validated, lift assets matched against specific movement requirements and 
schedules updated as necessary. For each movement requirement (force, equipment or 
supply) an approved operations plan designates lift assets, ports of embarkation, ports 
of debarkation, load dates and delivery dates. If the deployment plan cannot be 
supported, the decision makers must make one of the following decisions. They must 
allocate additional lift assets, adjust required delivery dates, delete low priority 
movement requirements or change the operations plan. 

2. Crisis Planning 

The second category is crisis planning. Crisis action deployment planning 
requires decisions on force commitment within an initial four hour period [Ref. 2]. The 
Supported Conimander must be able to develop a course of action which 1s acceptable 
in terms of his military objectives and is supportable by the transportation system 
[Ref. 3: p. 1-2]. The key question is, “Can the plan be supported?” Furthermore, 
guidance such as size of deploying forces, quantity of material, ports of embarkation 
and debarkation, required delivery dates or rates, and mode selection must be passed 
down to the Commands controlling land, air, and sealift assets. These Transportation 
Operating Agencies (MTMC, MAC, MSC) must initiate their planning processes in 


order to properly support the deployment in a timely manner. 
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For the analyst the key difference between deliberate and crisis planning is the 
time constraint for constructing a viable deployment plan. Any model that is used as 
an aid in the crisis decision making process must generate a solution in under four 
mours |Ref. 2]. 


D. INPUTS TO PROBLEM 
In order to coordinate land, air and sealift assets the JDA assembles data relating 
to deployability of all forces, equipment and supplies. 


e The first inputs are the movement requirements. A movement requirement 
consists of the following items: 


- Classification (passengers, equipment, fuel, etc.) 
- Weight and Quantity 

- Origin and Destination 

- Available Load Dates 

- Required Delivery Dates 


e The second inputs are a listing of the ports and their capacities. Capacity is the 
short tons per day of material that can be moved through the airport or 
seaport. 


e The third inputs are the lift assets. 
-  Tvpe and Quantity available 
- Capacity in short tons 
-  Utuihzation rate 
- Transit times 
Given these inputs the Joint Deployment Agency attempts to optimize strategic lift 


resources in support of deployments. 


E. RESEARCH TOPIC 

This study will develop a prototype model that can be used to allocate strategic 
sealift resources. The model will attempt to show that integer programming 1s a viable 
method for solving the sealift allocation problem. Furthermore, the model will 
incorporate intelligent reduction methods to restrict the number of equations and 
variables so as to reduce the size of the model. Without these reductions, the integer 


programming approach would not be computationally viable. 


Us 


It. BACKGROUND 


A. GEORGIA TECH RESEARCH GROUP 
The primarv research attempting to solve this large deployment problem has 
come from a group led by Professors John J. Jarvis and H. Donald Ratliff of the 


Georgia Institute of Technology. During the past five years their effort has been to: 


1. examine deployment planning in a crisis action environment from a 
modelling perspective; 

2. assess available methodology and modelling concepts for 

application to the crisis action environment; 

develop concepts and methodology for closure optimization; 

4. develop a system architecture within which these models would function 
[Ref 3p ae 


Lor) 


Jarvis and Ratliff describe four levels of models, which if developed, would 
function in a hierarchy. Decisions and assumptions made at the higher levels would 
guide and constrain decisions at the lower levels. Violations of these constraints could 
not occur unless the higher level modifies or changes the constraining decision or 
assumptions. The lower the level the more the detail involved in the planning process. 

The highest level is the closure planning level. The primary purpose of this level 
is to aid the decision maker in developing a general movement plan which will sausfy 
the military objectives and can be supported by the available transportation system. A 
general movement plan includes mode, port of embarkation (POE), port of debarkation 
(POD), assignment of movement requirements, timing of movements, degree of 
flexibility allowed at lower levels, and manner in which movement requirements can be 
split for transportation. The decisions made at this level are the most important 
because they guide and constrain all future decisions. [Ref. 3: pp. 7-10] 

The second level is the system loading, coordination level. Its purpose is to insure 
efficient utilization of the transportation system in carrying out the general movement 
plan developed in level one. At this level they search for and attempt to resolve 
problem areas, and develop more detail regarding movements. Additionally, it provides 
information and coordination necessary for transition from the top level to the detailed 


scheduling by transportation operating agencies in level three. [Ref. 3: pp. 9,11] 
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ie CoinGmlevelmMon ncrescerailed schedules for deployment are constructed. The 
three transportation operating agencies, MT MC, MAC and MSC are required to 
construct schedules to satisfy the specifications set forth by the first and second level. 
These specifications include movement requirements, suggested lift asset, POE, POD, 
and the required delivery dates. [Ref. 3: p. 9] 

Level four system is for monitoring the development and implementation of the 
deployment plan. This four level system is a dynamic planning system that provides for 


feedback, updates, and modifications as the plan proceeds. [Ref. 3: p. 10] 


B. GEORGIA TECH MODEL 

ie smcineenrusmmethe Georgia Pech researchers has been on level one, where 
the general movement plan is developed. They decided the best way to solve the 
deployment problem was to use decomposition. They broke the problem into two 
subproblems, a channel configuration and a movement requirement assignment 
problem. The problems are connected through a set of linking constraints. The 
decomposition method first generates the solution to the channel configuration model. 
With the link constraints fixed, the miovement requirement assignment problem is 
solved. The results of this model generate a linking constraint that 1s passed back to 
the channel configuration model which is solved again. This process is repeated until 
the solutions converge at optimality or it can be stopped at the user’s discretion if time 
is limited. [Ref 3: p. 45] 

In the search for appropriate solvers Jarvis and Ratliff have attempted to take 
advantage of the sparseness of the network generated by a deplovment problem. The 
movement requirement assignment problem can be formulated as a pure network 
Structure, therefore it can be best solved using a network solver [Ref. 3: p. 49]. For the 
channel configuration model they chose a solver for a network with side constraints 
[Ref. 3: p. 245]. Due to the large number of side constraints they may switch to a 
linear program solver in the future [Ref. 2]. The two problems are linked together with 


Benders decomposition method. 


C. MODEL PROBLEMS 

The current model is experiencing several problems. The solution procedure 
adopted for the model converges slowly and at times does not produce an optimal 
answer. Research is ongoing in an attempt to discover the source of the convergence 


problem. 
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When the procedure produces an optimal solution, it requires over six hours even 
for a medium size deployment problem. This 1s not acceptable for crisis planning. 
Current research is investigating a “hot restart” capability, aggregation of movement 
requirements, suboptimal stopping rules, a method to generate arcs as needed, and arc 
reduction methods. 

A third area of concern is the method of modelling sealift. The model assumes 
continuous flow variables even though cargos move discretely in reality. The concept 
can best be understood by likening the channel and its capacity to a pipe with a water 
passing through it at a given flow rate. 

The Georgia Tech research group makes a good argument for the channel 
concept and continuous flow rate when applied to airlift. The airlift cycle times are 
relatively small when compared to the time horizon and the delivery effect is 
“smoothed” over time. [Ref. 3: p. 35] However, the assumption of continuous flow 
variables for sealift 1s not realistic. Travel times are not small when compared to the 
time-step of the model. Furthermore, the assumption of a steady-state system cannot 
be made during the initial 20 to 30 days and it 1s questionable whether the steady-state 
system argument can be applied to each port in the later stages of a deployment plan 
(30-180 days). Finally, the simplifying assumption of continuous flow variables for 
sealift is not necessary. This study will show that ship departures and arrivals can be 


modelled discretely. 


D. OTHER RESEARCH 
In a related Masters thesis in Operations Research, CPT. Ko Sievens@aliew 
Naval Postgraduate School, Monterey, California constructed a model for determining 


optimal flow rates for air and sealift. [Ref. 4] 
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Hie rOpDEL DEVELOPMENT 


A. INTRODUCTION 

inet o4 Oe. Dantzig and DM. R. Fulkerson showed that the problem of 
determining the minimum number of tankers required to meet a fixed schedule could be 
modeled and solved as a linear network flow optimization problem [Ref. 5: p. 217]. 
The problem assumes that the load date, origin, and destination for each shipload are 
known, all ships have the same capacitv and speed, and all movement requirements are 
full shiploads. This study extends the Dantzig and Fulkerson tanker scheduling model 
to include time windows for loading. Instead of specifving a single load date, a time 
window from Earliest Load Date to Latest Load Date may be specified. The 
remainder of this chapter will compare the two formulations and present an example of 


the new one. 


B. MATHEMATICAL FORMULATIONS 
]. Dantzig-Fulkerson Model 


Indices: 
1 = Port of embarkation (POE) index represents the port of 
origin for a shipload of cargo 
] = Port of debarkation (POD) index represents the port of 
destination for a shipload of cargo 
a = Loading time period for cargo 
b = Unloading’ discharge time period for cargo 
Given Data: 
qT; = Loading time + transit time from POE 1 to POD | 
Wes = Unloading time + transit time from POD j} to POE 1 


1} 
Derived Data: 


ile § = Number of ships that load at POE 1 on day a 


Nbj = Number of ships that unload at POD j on day b 


Decision Variables: 


N aibj = Number of ships departing from POD j at time b to 
return to POE 1 at time a (reassignment of ship) 
Oe = Number of ships starting their schedule at time a from 
POEA 
Ypj = Number of ships completing their schedule at time b 
at POD } 
Formulation: 
Viaxinmize eZ 
Silsjecusl oc: 
ys Nail = Oe for all a,i (eqn 3.1) 
b.] 
Giger iy Ng for all b,j (eqn 3.2) 
a,i 
y Xie a ee (eqn 3.3) 
a1 al 
b4] b,] 
» Nibie 2 (eqn 3.5) 
a,i,b,] 


XXqi 20, Yj; eee 0 


aib} 
Constraints 3.1 require every shipload of cargo be moved. The cargo may be 


moved by a ship on the first leg in its schedule, XX,:, or a ship that is reassigned for 


al’ 
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an additional leg in its schedule, X,:,;. Constraints 3.2 insure all required shiploads of 


cargo are delivered. The delivery ei Demiie Jdst StOp Of a ship, Y pj: Ge One On several 
reassignments, Naibj within its schedule. Constraints 3.3 and 3.4 insure that there are 
sufficient number of ships to satisfy all pickup and delivery requirements. The 
objective function miaximizes the total ships reassigned, which is equivalent to 
minimizing the number of new ships used. 

Por aa Naibj to be feasible the arrival day a at POE 1 must be greater then or 
equal to the unloading day b at POD j added to the transit time between the two ports. 


merollows that if 


a<br ITi 


is infeasible and X Si@ed be mMRecero ZerOner, preferable deleted from 


Bes aid) aibj 
the problem. This condition is not an explicit constraint in the model, but rather the 


relationship is used to reduce the set of X,:1;’s to insure only feasible values are 


aib} 
defined in the model. 

The Dantzig and Fulkerson formulation is totally unimodular, therefore the 
optimal solution is an integral solution. 

2. Extended Formulation with Time Windows 

As stated earlier, the new formulation extends the Dantzig and Fulkerson 
tanker scheduling model to include time windows for loading. Instead of specifying a 
single load date, a time window from Earliest Load Date to Latest Load Date may be 
specified. This situation can be equivalently modelled using time windows for 
unloading or arrival dates. The time window would specify Earliest Arrival Date to 
Latest Arrival Date. With the addition of the time windows, the extended formulation 
is no longer totally unimodular, therefore you are not guaranteed an integral solution. 
To insure an integral solution it is necessary to use an integer program solver. The 
integer program solver enables the user to model ship flows discretely versus the 
Georgia Tech model which assumes continuous flow variables. Similar to the original 
formulation, there is a specified origin and destination for each shipload and all ships 
are considered equal in capacity and speed. 

Using one type ship in the model versus choosing between several different 
types of ship is acceptable at the strategic level. At this level the objective is to 


determine if a deployment plan is feasible. If the plan is feasible, a lower level will 
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implement the actual movement plan using more detail. At the strategic level cargo is 


broken into three categories; unit equipment, resupply and ammunition. For each 


category a specific ship type is designated. Unit equipment will make maximum use of 


roll-on roll-off ships. 


Resupply cargo and ammunition will use containerships and 


breakbulk cargo ships. The mix between the latter two is determined by the theatre of 


operations and the extent destination ports can handle containerships. The model can 


be run for each type ship and category of cargo. Any movement requirements not 


moved by the allocated lift assets in the first run, can be included with cargo in the 


next category during subsequent runs. 


Indices: 


b 
Given Data: 


TT; 


REQ aj 


Derived Data: 


POE, 


POD, 


Sal 


Port of embarkation (POE) index represents the pomrans 


origin for a shipload of cargo 


Port of debarkation (POD) index represents the port of 


destination for a shipload of cargo 
Loading time period for cargo 


Unloading/discharge time period for cargo 


Loading time + transit time from POE1 to POD j 
Unloading time + transit time from POD j to POE 1 


1, if shipload s is required to load at POE 1 on day a 
with destination POD | 


QO, otherwise 


1, if shipload s departs POE 1 


0, otherwise 


onl 


, if shipload s arrives POD | 


QO, otherwise 


-— 


, 1f shipload s loads at POE 1 on day a 


QO, otherwise 
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Dsp; = 1, if shipload s unloads at POD j on day b 


OQ, otherwise 


LINK gi = 1, if shipload s goes from POE ito POD ; 


0, otherwise 
Decision Variables: 


N aibj = Number of ships departing from POD } at time b to 


return to POE 1 at time a (reassignment of ship) 


LO = Number of ships starting their schedule at time a from 
POEW 

Y bj = Number of ships completing their schedule at time b 
at POD | 

eae = 1, if shipload s departs POE 1 on day a 


QO, otherwise 
Formulation: 


Maximize Z 


puoject To: 
i) en Oe sai for all a,i (eqn 3.6) 
b,] S 
eink bye oe sai for all b,j (eqn 3.7) 
al Sia 
Dav = | for all s (eqn 3.8) 
al 
DeXXa + Z =)! POE. (eqn 3.9) 
el S,1 
» Mp ae ~ POD, (eqn 3.10) 
b,j S,] 
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ex 
aliicloy 


aibj — =Z Coals) 


KX, 20) Vy, Scere 


Maes Ovorel 


Bay) 


Constraints 3.6 require everv shipload of cargo be moved. The cargo may be 


moved by a ship on the first leg in its schedule, XX,:, or a ship that is reassigned for 


aye 


an additional leg in its schedule, X Constraints 3.7 insure all required shiploads of 


aby’ 
cargo are delivered. The delivery may be the last stop ofa ship, Y pj or one of several 


reassignments, X within its schedule. Constraints 3.8 are generalized upper bounds 


aibj 
that pick one load _ from the time window of available load dates. Constraints 3.9 
and 3.10 insure that there are sufficient number of ships to satisfy all pickup and 
delivery requirements. The objective function maximizes the total ships reassigned. 
the movement requirements. 

Oftentimes, in a crisis deployment, the objective function is to minimize the 
time required to complete all the movement requirements with the available resource. 
In this model, the movement requirements have specific times they must depart in 
order to arrive at the destination and satisfv the Commander’s military objectives. 
Therfore, the optimal solution of this model determines the minimum number of ships 
required to satisfy all the movement requirements. Given sufficient assets the model 
will determine if there is a way to satisfy the requirements with fewer assets. If there 
are not sufficient sealift assets available, the Commander must allocate additional lift 
assets, adjust required delivery dates or change the operations plan. 


As in the Dantzig and Fulkerson formulation, for an X to be feasible, 


aibj 


Figures 3.1 and 3.2 are network diagrams which depict the variables used in 
the formulation. The general network diagram with all possible arcs is shown in Figure 
3.1. Figure 3.2 illustrates the effect of variable reduction methods that will be 
discussed in Chapter IV. The problem that generates these network diagrams 1s 


discussed in Section C. 


m, \ om, / 
POEs 
oS 





Eng OO 
Fe O O\O\O We 


TIME PERIODS 


figure 3.2) Network diagram with reduction methods incorporated. 


The new ships, XX,;, are arcs with tail at the source and head at node a,i. The ships 


that stop, YVibp are arcs With tail at node by and head at the sink. The reassigned ships, 
Said) saas 
be represcnied as an arc with tail at origin node a,i and head connected to the correct 
shy L. 

The kev difference between this formulation and the original Dantzig and 


Sorceress Witimtaieat node Dypand head at node at. Phe binary variable M 
desunation, node b,j, where D 


Fulkerson formulation is the introduction of the binary variable Myo). Por example, if 


shipload [is available to load at New York on day 1, 2 or 3 there will be an My; for 
eee dticeemand }) — NN) (See fieure 3.2). Constraints 3.8 select the one 


i 


Constraints 3.6 insure that for every arc M..; leaving node ai there is an arriving arc 


Naibj 


Ciitencprescitsmelespest loud dame irony tlie time window of avarlable dates. 


or arc NX, as shown in igure 3.3. 






oe. 
vA 


PST 3 me OMsURMIts 3-0. 


Constraints 3.7 insure that for every arc My...; arriving at node by there is a departing 
ee \aibj 
Mite bingiy variable Vi 


Or are Yj as shown in Figure 3.4. 
Seomirectwiitned tO «2 nove my, therclore it 1s 


Sal sb} I. The 
merece node ty ismlcternumed by carefully delining a subset for cach index $s.) and a. 


Sal 


necessary to find a method to link M eomuileecorrece Node, wiiere 1 


The righthand side of Constraints 3.7 is now rewritten more preciscly as follows; 


S Ss voM 


pos ee a St 


=e Wey Glee > Z 
SEDs, 1EPOER, aca 


lor all b,j. 


tJ 
Lee 





Figure 3.4 Constraints 3.7. 


The above notation, using GAMS inspired shorthand, reads, “Summing over all s such 
that Dobj > 0; Summing over all t such that POE. > 0; Summing over a such that a 


= b- Tye Examples of Doty and POL: are found in Section C2. 


C. EXAMPLE PROBLENI 
To illustrate the model. consider the following numerical example for Which the 
corresponding network diasran issdepictedstiim (etnies wet (ienemse 
2 POLS ie as 
PODs,) = Eng, Ir 


e 3 shiploads 


tr 


© 10 tintetpeneds: 
1. Given Data 

There are three tables of data required to run the model. [In the cxample we 
assume the transit times are equal in both directions, qj; = TT. Therefore, the model 
input will consist of two sets of data. Table 1 lists the time it takes a ship to load or 
unload its cargo and travel from port 1 to port j. The travel times are purposely kept 
small. 

The Movement Requirements Table (REQ) lists cach shipload. port of 
embarkation, available to load dates, and port of debarkation. Vollowing the notation 
used in the GAMS modeling language [Ref. 6.7,8.9,10] Table 2 reads, “Shipload I is 


available for loading in New York on days | through 3 with a destination of England.” 


TA 1 


Vij ives (DAYS) 


Eng ene 
NY, l 2 
Jax 3 2 


Each shipload must have one POE i and one POD j. Our GAMS implementation of 
the model contains a check for this condition. If an error is detected the program will 


stop immediately and print several error messages ta help locate the nustake. 


TABLE 2 
MOVEMENT REQUIREMENTS 


Eng ee 
SOlENeY TOT" 103 l 
S02. sv. 103" 104 
503.Jax.T04* 105 





2. Derived Data 

After checking for errors in the REQ Table several other tables of information 
are derived from the REQ Table. The information 1s derived through a process that 
Kendrick and Meeraus call projection. The GAMS modeling language stores data as a 
relational database. “Projections can be used to answer queries by projecting the n 
dimensional data in a relationship onto an n-1 dimensional or smaller space.” 
[Ref. 8: ch. 19, p. 3] The derived data is used to reduce the size of the model and 1s 
used in the procedure that determines the schedule for each ship. Listed in the 
following tables are examples of the derived data. 

Table 3 is an example of projecting the four dimensional data in the REQsiaj 
Table onto the three dimensional space of the Oe Table. Table 3 reads, “Shipload #1 
is available for loading in New York on days 1, 2 or 3. Shipload #2 is available for 


loading in New York on days 3 or 4.” 


oes 


TABLE 3 


Example of O..; 


NY Jax 

S01.TOL l 
S01.TO? l 
S01.T03 l 
| S02.T03 l 

S02.T04 l 
S03.T04 l 
S03.T05 l 


Table 4 gives possible dates for unloading. For example, the first five rows 
read, “Shipload #1 unloads in England on day 2, 3 or 4. Shipload #2 unloads in 


Prance on dav 01.6. 


TABLE 4 
Example om Dev 


Eng Fr 
SUE i2 | 
SO1.T03 | 
S01.T04 l 
SO. OS | 
$02.T06 | 
S05.007 | 
Soe. 008 | 


Table 5 determines the connection between the origin and destination. It 
reads, “Shipload #1 departs New York with a destination of England. Shipload #2 


departs New York with a destination of France.” 


TABLE 5 
leqesieajelle ert LINK gi 


Eng lee 
SUlEN i 
SO2IN Y ] 
S03.Jax J 


Table 6 reads, “Shipload #1 and #2 load in New York and Shipload #3 loads 
in Jacksonville.” 


TABLE 6 
Example of POE,, 


oa Jax 
S01 ] 
S02 J 
S03 I 


Table 7 reads, “Shipload #1 and #3 unload in England and Shipload #2 
unloads in France.” 


WADIEE 
Example of POD, 


Eng eye 
S01 J 
S02 ] 
S03 ] 


aT 


IV. MODEL REDUCTION 


As formulated above, the problem includes a large number of variables and 
constraints even for a medium size problem. This, in practice, means that to obtain a 
solution would require a large amount of computer memory and cpu time. 
Additionally, for a large problem, the model may be too large for most computers. 
Therefore, it 1s important to develop intelligent methods to reduce the number of 
variabies and equations used in the model. This chapter will discuss several reduction 
methods that make the new model computationally viable. 

The General Algebraic Modelling System (GAMS), developed by Alexander 
Meeraus and Tony Brooke, provides modelers with a tool to eliminate all variables 
unnecessary for the problem in an efficient manner, which reduces the size of the 
resulting optimization problem. An important feature is the S$ operator, which 
translates to “such that’. The §$ operator can be used to restrict the Gomameson 
variables, reduce the number of equations generated, or specify that the summation 
should be only over a defined subset of indices. The remainder of this section will use 
GAMS notation to explain model reduction methods. 

Constraint 3.8 provides a excellent example to illustrate the use of GAMS 
notation. Constraint 3.8 is referred to as LOAD(s) in the GAMS formulation. LOAD 
is the “name” of the equation and the set of shiploads s is the “domain” of this 


equation. If Constraint 3.8 were written for all s it would appear as follows; 
LOAD(s) ... SUM((a,1), M(s,a,1)) =E= I. 


This reads, “Generate a LOAD equation for each s. In each equation the sum over a,1 


for all ee 


feasible set for the binary variable M,,;. The smaller the set of binary variables, the 


equals one.” In the formulation it is important to accuratelv define the 
faster the problem can be solved. Feasible M,,.’s are possible only when there is a 
shipload s available for loading at POE i on day a. Therefore, the set for defined 
variables M,..; 1s restricted by using the data derived for the O,,; Table. Constraint 
3.8 is nOW Written; 


LOAD(s) ... SUM((a,i), M(s,a.i) § O(s,a,i)) =E= 1. 
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sins eiedds; sum over a, all M_.: 


, such that s,a,1 belongs to the set defined by 
Oca) > 0. 


A. VARIABLE REDUCTION 
In 1954 the Dantzig-Fulkerson formulation was designed to be solved manually 


using the simplex algorithm. They identified the feasible subset of X by creating a 


aibj 


tableau and crossing out all cells containing infeasible X With GAMS notation 


aibj’ 
this can be accomplished with the following lines; 


SE TerPeASIB IEE feasible reassignments for X; 
FEASIBLE(a,1,b,j) = YES S (ord(a) - ord(b) GE T(i,j)). 


This reads, “A subset {called FEASIBLE) of the indices a.i,b,j 1s defined for only those 
combinations of a,i,b,j such that the ordinal position of a minus the ordinal position of 
b is greater than or equal to qi; {a-b 2 Tipe With this set defined it is now possible 
to generate equations that include only the feasible variables. For example, Constraint 


3.11 becomes 
SUM(a,1,5,}), X(a,1,b,j) S FEASIBLE(a,i,b,j)) = Z. 


This reads, “Sum over a,i,b,j all Naibj SDenmtyatwa lo. isameene subset FEASIBLE ... ~ 


A second method to reduce the set of X aibj 1S 


SET PORT 
PORT(a,i,b,j) = YES S SUM(s, O(s,a,i)) S SUM(s, D(s,b,j)). 


Bach X 


mun) N. 


aibj is an arc with tail at node b,j and head at node a,1. Set PORT forces every 


aibj b° have a tail where a ship has delivered cargo. If no ship had delivered 
cargo to that port there would be no ship available for a return trip. Furthermore, Set 


moo icrces every arc X to have a head where there is a shipload ready to be 


aibj 
loaded. There is no reason to send a ship to a POE if it 1s not needed at that location. 


A third method to reduce the set X is based on common sense and 


‘ aibj 
experience, and is controlled by the user of the model. 


SET RESTRICT 
RESTRICT(a,i,b,j) = YES S (ord(a) - ord(b) LE (T(i,j) + SLACK))). 
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The value for variable SLACK is set by the user of the modell) A ship mustuge 
reassigned to a POE 1 within a specified length of time (transit time + slack time). If 
the ship is not reassigned within the time period it will stop. This is indicated by a 


value for Yy;. A ship will always be reassigned if possible because the objective 


function ee reassignments. A requirement outside the time period will use a 
new ship. When looking at ship schedules it will be obvious if two schedules could be 
satisfied by one ship. The final time period in one schedule will be less than the initial 
time period in another schedule. The use of set Restrict 1s acceptable because onemeam 
assume a ship returned to its homeport if not reassigned in a timely manner. In the 
future the ship would be available at a new location when needed. 


In the formulation all three set reduction methods are combined into 


SET aOR X values allowed in the model 
OK(a,1,b,)) = YES $§ FEASIBLE(a,i,b,j) 
S PORT(a,i,b,j) 
SeRESTRICT({at6,)). 


Constraint 3.11 becomes 


SUM((a.i.b,j), X(a.i,b,j) S OK(a,i.b,j)) = Z. 


B. EQUATION REDUCTION 

In addition to reducing variables it is important to reduce the number of 
equations in a large model. Constraint 3.6 is specified above for every possible 
combination of a and i. However, it is only necessary to generate equations when 
there 1s the possibility of some shipload s loading at POE 1 on dav a. With GAMS, 
we can reduce the number of equations by generating Constraint 3.1! (or POECAP(a.1) 
as it 1s called in the GAMS formulation) only when needed. We do this by applying a 


S operator to the equation definition; 


POECAP(a,i) $ SUM(s, O(s,a,i)) « 
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This is interpreted , “Generate an equation for each (a,1) pair such that the sum over s 
of O,,- is non-zero.” This condition will be true and hence the equation will be 
generated only when at least one shipload is eligible to leave POE 1 on day a. The 
remainder of the equation can be seen in Appendix A. 

Similarly, Constraint 3.7 should onlv be generated when there 1s the possibility of 
a shipload s unloading at POD j on day b. Therefore Constraint 3.7 , which is referred 
fopas | ODCAP(0.};) in the GAMS formulation, is restricted as follows; 


PODCAP(b,j) $ SUM(s, D(s,b,j)) «.. 


The interpretation 1s similar to Constraint 3.6. 


C. GAMS SOLUTION TO THE EXAMPLE PROBLEM 
The optimal solution of this formulation consists of four sets of variables. Also 
listed is the solution of the example problem. The network diagram for this optimal 


solution is shown in Figure 4.1. 


oe XX. 


i 4 schedule is constructed for each new ship. 


XA(TOLNY) = 1 
XX(T04.Jax) = 1 


e. Naibj is a reassignment of a ship from POD j dav b to POE 1 day a. 


Me iGeeny.f027Eng) = 1 

e, Y pj a schedule 1s completed for the ship at POD j on day b. 
SACO IE se)) 
Y(T07.Eng) = 1 

«. M,,; 1s the connection between POE 1 and POD }. 


M(SO1.TOL.NY) = 1 
M(S02.T03.NY) = 1 
M(S03.T04.Jax) = 1 


The output is printed to a file which is then input to a Fortran program, which 


constructs a schedule for each ship. 
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Iigure 4.1 Optimal solution to the example problem. 


V. SCHEDULE PROCEDURE 


Dantzig and Fulkerson developed a procedure to turn their solution into a 
schedule for each ship [Ref. 5: pp. 219-220]. We present a new procedure, implemented 
in a FORTRAN program called SCHED, that is an adaptation to accommodate the 


binary variables M Given integral solutions to the four sets of variables from the 


Sal’ 
GAMS model, SCHED will construct a schedule for cach ship. SCHED has two 
major components. Part | reads in the data from an output file generated by the 
GAMS model. If the solution is non optimal, infeasible, or unbounded the program 
will stop and print an error message. If the solution is optimal it goes to Part 2 which 


contains the following procedure for determining a schedule for each ship. 


Ship Schedule Procedure 
Input: Integer Values for Xaibj? XXai, Vp}: and olga 
Output: Ship Schedules 


For all xx,,; = 0 
feoclect any a.) 


Erm Depart 1, Day’, a) 


Select any s such that my, = | 
OR ete a 
Meme a > 


Select anv j° such that LINK. 4 = | 
b' = Ty +a 
Print (Arrive’, }’, “Day”, b’) 
Do 
{X = { Xaib’y’ for all a,1: Xaibj’ Zz (0) ) 
While (X + 9 ) 
Select any a’,1 such that Xaiby’ © x 
Pimimieecoinen), Day ab) Arveem, Day’, a) 


Select any s such that m.:,- = | 


ac oil = a Toy Ie = ] 
Mo’ 4’j’ = IN’ a 7’ - l 


Select anv }’ such that LINK Gy’ = | 


b = ey a 
Print Atiive ice bay >) 


Print (“stop”) 


Ybj’ = a - l 


ne 


It is possible to generate several different schedules from a GAMS solution. 
However, the total number of ships utilized is not changed. The automation of the 
process, of converting GAMS output to SCHED input, makes this ship schedule 


procedure very easy to use. 


A. SCHEDULE OUTPUT 
The output from SCHED 1s a listing of the schedules for each ship. The output 


from the example problem 1s shown in Table 8. 


TABLE §& 
OPTIMAL SHIP SCHEDULES FOR EXAMPLE PROEMI=™ 


Ship #1 POR PERTOD SHIPLOAD 
Depart Ne l ] 
Arrive Eng Z 
Depart Eng ps 
Arrive NY 3 
Depart NY 3 2 
Arrive eae 5 
stop 
Ship #2 PORT PERIOD SHIPLOAD 
Depart Jax 4 3 
Arrive Eng 7 
Stop. 
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B. TIME HORIZON 

A method that could be used to reduce the size of models with long time 
horizons is to increase the size of the time period. This reformulation would require 
the interpretation of the data, indices and variables to all change. Changing the time 
unit from 1 day to 2 or 3 days for large sealift model would be reasonable. However, a 
time unit of four days or more is difficult to justify. The longer time periods give less 
resolution to the model. Furthermore, the solutions are more general and harder to 


interpret. 


VI. RESULTS 


A. IMPLEMENTATION 

GAMS was initially used to develop the model. During the development stage it 
saved time because it enables the user to input models using algebraic notation. This 
permits one to make changes and experiment with different concepts without becoming 
bogged down in programming a model generator. The error messages generated from 
mistakes or errors in logic were helpful in debugging and correcting problems. 
Additionally, the use of the $ operator, translated to “such that”, enabled the model to 
restrict the domain of variables, reduce unnecessary equations, specify subsets for 
indices and perform data transformations. 

Once the problem was formulated it was necessarv to choose an integer solver. 
GAMS is able to interface with the integer solver ZOOM, therefore GAMS was also 
used to generate the model. 

ZOOM. XMP is a mainframe and PC computer code for solving zero/one mixed- 
integer programming problems developed by Roy E. Marsten at the University of 
Arizona. [Ref. 11] ZOOM solves the model generated by GAMS and its output 1s 
written to a formated file which serves as input to the Fortran program, SCHED. 

SCHED, the Fortran program that constructs schedules for each ship is written 
in FORTRAN 77 and compiled by the IBM VS FORTRAN complier. It takes input 
from GAMS'‘ZOOM and constructs a listing of schedules. 


B. COMPUTATIONAL RESULTS 

Modeling the allocation of strategic sealift using integer programming in 
conjunction with intelligent reduction methods is viable. Integer programmung enables 
the user to model discrete shipments versus modelling a continuous “flow” of cargo 
into sea ports. Problems involving 120 binary variables have been solved with 
GAMS/ ZOOM in 3.6 minutes on an IBM AT with math coprocessor running at 8 
megahertz. 

Intelligent reduction methods are critical to solving large problems. Restricting 
the domain of variables and deleting unnecessary nodes and arcs in the model reduced 
the size of the problem. The following two examples illustrate the impact of the 


reduction methods. 


Example Problem (Mata given in Chapter 3C): 
me FOES 
ee 2 POWs 
@ 3 shiploads 
e 10 time periods 
e SLACK = 5 days 


Before After Percent 


Reduction Reduction Reduction 


Constraints 46 19 9 
Continuous Variables 440 18 .96 
Binary Variables 60 7 88 
Non-zero Elements 1260 62 a) 


Example 2 (Data given in Appendix B): 
wee POEs 
e 5 PODs 
e 30 shiploads 
e 60 time periods 
@ SLACK = 5 days 


Before After Percent 


Reduction Reduction Reduction 


Constraints 633 203 .68 
Continuous Variables 90,600 589 99 
Binary Variables 9000 96 99 
Non-zero Elements 298,000 1885 99 


This significant reduction in the number of variables, equations and nonzeroes 
makes the difference between success and failure of the integer programming 
formulation. As the time horizon is extended the user controlled restriction of Naibj > 
via the SLACK parameter has an even greater impact. The SLACK parameter 
requires that a ship be reassigned within a specified length of time. If it is not 


reassigned within the time interval it will stop. Too small a value may over restrict the 


problem and result in poor solutions. Too large a value will not effect the solutions, 
but will fail to reduce the size of the problem. Experience and careful consideration 


must be used in choosing an appropriate value for SLACK. 
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Vile CONEGEUSTIONS 


The transportation community has expended a tremendous amount of time and 
effort attempting to efficiently allocate lift assets. In the past many of the efforts have 
been disorganized and geared to specific problems. The Unified Transportation 
Command, which will soon replace the Joint Deployment Agency, has a goal to insure 
that the data bases and computer models will be coordinated at every level [Ref. 1: p. 
44]. Under the guidance of the Transportation Command the research effort in large 
deployments will continue. 

Results of this study indicate that integer programming solvers used in 
conjunction with intelligent model reduction methods can find solutions to deployment 
problems. These solutions can be used to aid in the planning and execution of 
deployments. Furthermore, integer programming with the resulting discrete ship flows 
is a Viable alternative to the Georgia Tech method of modeling sealift as continuous 
flow variables. 

There are several areas for future research that can improve the model. 


e Generate sensitivity analysis information on the time windows. The insights 
gained would be invaluable both in refining existing deployment plans and 
developing new deployment plans. 


¢ Develop the capability to generate alternate optimal solutions. Alternate 
solutions would allow the deployment planning process to be more flexible. 


e Develop a model that allows a mix of ships with different capacities and speeds. 
The model would then determine the best ship for a specific movement 
fe quirement, 


e Pernut multiple loading (POE) and discharge ports (POD) for a single shipload. 
A ship could pick up cargo at two or more locations and deliver to several 
discharge ports. 


e Rather then aggregate movement requirements into shiploads, it might be better 
to match individual movement requirements against an asset. This idea enables 
the model to consider cargo compatibility and realistic groupings of movement 
requirements. 


A potentially promising approach that can be used to implement many of these 
ideas 1s to use column generation [Ref. 12]. A ship schedule generator could portray 
realistic shiploads and include any restrictions on ship type, capacity, speed, port access 


Or cargo compatibility. In addition to answering the strategic allocation problem, a 


column generation formulation would benefit the Military Sealift Command (MSC) in 
the scheduling of individual ships. The increased realism and the numerous parameters 
and constraints that can be modelled with an intelligent generator would increase user 


“faith” in the solution. 
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APPENDIX A 
GAMS FORMULATION 


offupper 
offsymlist offuellist offuelxref offsymxref 
sontext 


Sof ur EXAMPLE PROBLEM 


Remerence: 

1. “Minimizing the Number of Tankers to Meet a Fixed Schedule", 
G.b. Dantzig and D.R. Fulkerson, Naval Research Logistics 
PucdacomemeVvOlweieetrs, Pp.2l/=222, September, 1954. 


Comment: This 1s a prototype model which can be used to allocate 
strategic sealift resources. It is a GAMS-Based Integer Program that 
extends a classic network flow optimization model developed by Dantzig 
and Fulkerson. Using intelligent reduction methods, the formulation 
reduces the number of constraints by 60-70% and the number of 
variables and nonzero elements in the matrix by 90-99%. 


Sofftext 


SETS i POE /NY,Jax/ 
5 POD maser 
Ss SHIPLOADS/S01%*S03/ 
TI TIME INTERVAL/TO1*T10/ 


eras (a,b,1I); 


TABLE T(i,3) TRAVEL TIME BETWEEN PORTS 
ENG FR 
NY it Z 
Jax 3 2 
TABLE REQ(S, 2.ana) REQUIREMENTS TABLE 
ENG FR 
S01.NY.TO1*TO3 ] 
S02. NY. fO0s7104 il 
S03.Jax.TO4*T05 1 
Sontext 


Hie spabamecen and scalar Valticswbelow are used to insure no 
mistakes were made when inputting the data into the REQ Table. 
1f a mistake is detected an error message will be written. 


a Srecee = 1, implies 1 destination was given (GOOD). 
Pymecheci smo 2, lmpreles PZ destinations were even (ERROR). 
€) if check(s) for a shipload is not listed t 

destination was given (ERROR). 


Sofftext 


PARAMETER CHECK(s) insure each Se ae unique destination; 
CHECK(s) = SUM(j $ (SUM((i,a),REQ(s,i,a,j))).1); 


SCALAR ERROR shiploads do not have a unique destination; 
ERROR = 1 § (SUM(s $ (CHECK(s) NE 1),1)); 


DISPLAY $ (error) "One or more shiploads does not have a unique POD."; 
DISPLAY $ (error) CHECK; 

ABORT $ (error) "Find the error and correct the REQ Table."; 

Sontext 


The below parameters are generated from the above REQ Table. 
They are used to restrict the size of the model and are used 


at implies no 


4] 


in the algorithm that determines the schedules. 


sofftext 

PARAMETER 0O(s,a,i) shipload a loads at 1 at time a; 
O(s,a,;1) = 1 S$ SUM@WREO sae 

PARAMETER  D(s,b,j) shipload s unloads at j at time b; 


LOOP Gone Ccy >) I= 
1S SUMGREOCS Rio TG ean oe 


PARAMETER  LINK(s,i,j3) | shipload s goes from i to j; 
LINK(s,1,3) = 1 S SUMCaRPHOCs i174 joe 


PARAMETER POR Cs)1) shipload s loads at POE 1; 
POE (s,1) = 1 $ SUM(a,0(s,a,1)); 


PARAMETER POD(s,j) poop eee s unloads at POD j; 
POD (5,9) =) 1a Gest cae 


Sontext 


The below sets reduce the size of the problem by deleting 
infeasible and unwanted variables. 


SOLrmatexst 

SET FEAS@ELE feasible reassignments; a 
FEASIBLE(a,i,b,j}) = YES $ (ord(a) - ord(b) GE T(i,j)); 

SET EORT include only necessary ports; 


PORT(a,1,5,]) = YES SUM(s,0(s,a,1)) . 
S$ SUM(s,D(s,b,j)); 


SET SREALI SUG Gonnections for specutresemie, - . 
REALISTIC(a,1,b,j3) = YES $ (Cord(b) = ord ajmmemte. jer 


Sontext 


The value for SLACK is set by the user of this model. A ship 
Will be reassigned to a port within a specified length of 
time(travel time + slack eie) if ete sont cannes wes 
reassigned wt tmethe Speaenee Time period the sitpewiecece. 
a4 requirement outside thé time periodhwi luce sane veo 

When determining ship schedules tt will bemebyvicus ieee 
schedules could be combined inte onegsehedtle. Yampa 

a schedule will be less than XX(a,i) in another schedule. 


Sofftext 
SCALAR SLACK sillack time for Ship 5/- 


SE SRE Sees user can reduce the number of reassignments; 
RESTRICT(a,i1,0,3) = YES S$ (ord(a) = ord(b) LE 92,7) eciaee 
SET OK X values allowed in model; 


OK(a-3.5.7) = VEouo EAs IBLE (sone) aan 
$ RESTRICT(a,i,b,j) 
S$ PORT(a,1,b,}3 


e 
f 


VARIABLES ss | . | | 
X(a,1,5,))  # Reassignments from ) ateeime 5 Coe meeeamees 
hota =) ? Ships starting cd eselmo was: soles! 

Wb) # Ships stopping at time b at j 

Z total return ate 

Mi say) binary variable to choose best load date 
POSITIVE VaRiApLise es. 75-. 

BINARY VARIABLE, 

aioe oe = 1; 

*M.lo(s,a,i) = 0; 

EQUATIONS 

Sas POE suppl 

PODCAP(b, j POD deman . 
START Eotal ships starting 
FINISH total ships finishing 
LOAD(s) GUB for mae 

ANS Objective function; 


POECAP(a,i) S$ (SUM(s,0(s,a,i))).. 
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SUM((b,j), (a,i,b,5) $ OK(a,d,b,5)) + ¥K(a/4) 


Ulm Comil(sraet) 6 Ols,4,1))- 


p P(b S D ‘ ee 
Per 29 Nee sos HS ay, ee 3) 


Se 3 POEs | 
SUM(i S$ POE(s,1), 
SUM(a $ REALiSTIC(a,i,b,j), 
M(s,a,1)))); 


meemb(s) .. SUM((a,1),M(s,a,i) $ O(s,a,i)) =E= 1; 
START .. SUM((a,i) $ SUM(s,O(s,a,i)), XX(a, i)) ¢ 
E= sil (s, deer OF ( San); 
m@ifon.. SUM((b,3) S$ SUM(s,D(s,b,3)),¥(b,j)) + 
=E= dunt ts, Wie EOD Gs ya)) > 


Meee UM (ant, )4(a,1,5,j3) S OK(a,1,b6,3)) =E= Z; 


MODEL DFMOD/ALL/; 


OPTION solprint=off,sysout=off; 
OPTION limrow=0, limcol= 0; 


Sontext 


These options format the output for a fortran program that 
produces ship schedules. 


Sofftext 


SETiLON X:3: 
BEIION XX:3: 


Oe: 
GC: 
BeetON Y:3:0: 
Oe 
CO: 


SETION M:3: 
SETLON LINK: 
ee Tton T:3: 


SE 


sl: 
1: 
1; 
sul 
ale: 


~e O-- ~e bs. 


XSOLVE DFMOD USING LP MAXIMIZING Z; 
BeevE DEMOD USING MIP MAXIMIZING Z; 


Sontext 


These values are used to set the limits of DO Loops within the 
BORBeban program, SCHED. 


Sofftext 

SCALARS XCNT number of X values 
XXCNT number of XX values 

Vent number of Y values 

MCNT nunper sor Mey cmiles 
are numper of possible connections; 

XCNT = SUM((a ome Lia, 1, bg) gt); 

KXCNT = SUM me We peeide, 1) 1) 

YCNT = SUM b.4) ae L(b, i) 1}; 

MCNT = SUM((s,a,i) $M ie Boe a. 

TCNT = SUM((i,3$,1 


Dior lay XCNT, XXCNT, YCNT , MCNT, TCNT; 
Petula ket, D7 LINK, A.L,Y.L,M.L; 


APPENDIX B 
INPUT FOR EXAMPLE 2 


Der uereive 
5 POES, 1 = NY, Chieleeesca. Si 
5 PODs, | = Eng, Gergrwanor ap 
30 shiploads 


60 time periods. 


Table 9 and Table 10 are the input data for Example 2 


dm 


TAREE 
qj (DAYS) 
Eng Ger Nc Kor 
NY 6 6 6 i 
Chr i d) tl 18 
Jax 8 8 8 18 
Sea 20 20 20 8 
Se 19 19 19 2 


-4 


Jap 
18 
17 
17 


SUNY. TOL" 103 
See NY.1O01* 103 
SUs NY .106* 10s 
Sd NY.119* 12) 
SUSEN Yel 257127 
S06. yn 1367140 
Umea) 50-153 
S08.Chr. TO02*T04 
SOe-ehr. Fos 110 
so.Ghr. 109" 111 
Sle C@hr.'24*T 26 
ple Cnr. 1267 125 
$13.Chr.T44*T46 
$14.Jax.T02*T04 
Sle oax027 105 
$16.Jax.T12*T14 
Sy iax23* 7 25 
Slo wax. 397 142 
$19.Jax.T46*T48 
Seiesiel03* 105 

SZ US Te 

See S24 126 

SESE Sed Bead Bers: 

Sees 49" 15) 

Syersca, 101103 
S20 sea Ol £03 
Sz eoca. 110° 112 
ec occas l2- tld 
ee sca. | 25" 127 
$30.Sea.T43*T45 


TABLE 10 
MOVEMENT REQUIREMENTS 


Fr 


Kor 


Jap 


ad 


Can 


~ 


le 
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